Distributed Control of Angle-constrained Circular Formations 
using Bearing-only Measurements 

Shiyu Zhao, Feng Lin, Kemao Peng, Ben M. Chen and Tong H. Lee 



Abstract 

This paper studies distributed formation control of multiple agents in the plane using bearing-only measurements. It is assumed 
that each agent only measures the local bearings of their neighbor agents. The target formation considered in this paper is a circular 
formation, where each agent has exactly two neighbors. In the target formation, the angle subtended at each agent by their two 
neighbors is specified. We propose a distributed control law that stabilizes angle-constrained target formations merely using local 
bearing measurements. The stability of the target formation is analyzed based on Lyapunov approaches. We present a unified 
proof to show that our control law not only can ensure local exponential stability but also can give local finite-time stability. The 
exponential or finite-time stability can be easily switched by tuning a parameter in the control law. 

Index Terms 

Bearing-only measurement, Circular formation. Distributed control. Finite-time stability, Lyapunov approach. 



I. Introduction 

Distributed formation control of multiple agents has been investigated extensively in various settings. We here would like 
to highlight two aspects that are important to characterize a formation control problem. 

The first aspect is what kind of information each agent can obtain from their neighbors. In order to realize distributed 
position control of multiple agents, it is commonly assumed that each agent can obtain the (global or relative) positions of 
their neighbor agents through wireless communication. It is interesting to note that the position information consists of two 
kinds of partial information: distance and bearing. Formation control merely using partial information has become an active 
research area in recent years. The work in |[T], Q addresses formation coordination of mobile agents when each agent can 
only measure the distances to their neighbors. Formation control using bearing-only measurements has been studied in Q-Q. 

The second aspect is how the target formation is defined. Conventionally target formations are defined by specifying global 
or relative positions of agents. It is noticed that angles and inter-agent distances can also be used to define a target formation. 
The term angle as used here refers to the angle subtended at one agent by its two neighbors. In recent years, control of 
distance-constrained formations has received much attention p0)-p6). Control of angle-constrained formations has also attract 
some interest very recently ||3)-||7), | [T7| . For distance-constrained target formations, any rigid body transformation (i.e., rotation 
and translation) over the entire formation will not change the inter-agent distances. If the target formation is defined by angles, 
in addition to rigid body transformation, the target formation will also be invariant to scaling. Furthermore, without parallel 
rigidity constraints ||7), p7) , |18|, the angle-constrained target formation will not be affected by any edge parallel motion 
either. 
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We now characterize the formation control problem studied in this paper from the above two aspects. Our work considers 
formation control using bearing-only measurements. That is motivated by an important type of sensor: camera, which inherently 
is a bearing-only sensor and has been widely applied in many control-related tasks. Vision-based formation control |[8), Q, 
[ [T9] is a potential application of our research work. In this paper, it is assumed that each agent can only measure the bearings 
of their neighbors in their local coordinate frame. It should be noted that the bearing measurements are not used to estimate 
any agent's position. The control is implemented directly based on bearing measurements. As a consequence, we can expect 
that the control is only able to handle angle-constrained target formations. Specifically our work considers angle-constrained 
circular target formations, where each agent has exactly two neighbors. In the target formation, the angle subtended at each 
agent by their two neighbors is specified. At last, our work makes no parallel rigidity assumptions ||7), p7) , p8) of the target 
formation. 

We propose a distributed nonlinear control law that stabiUzes angle-constrained target formations merely using local bearing 
measurements. The proposed control law is inspired by the work in ||3), Q, Q which addressed distributed control of triangular 
and quadrilateral formations. An attractive feature of the control law in ||3), Q is that the global stability can be proved based 
on the Poincare-Bendixson theorem. In contrast to previous work, we consider circular target formations with an arbitrary 
number of agents and the stability by the proposed control law will be analyzed based on Lyapunov approaches. Our work 
also involves finite-time formation control which has received some attention recently pO[ , pT) . Besides fast convergence, 
finite-time control can also bring benefits such as disturbance rejection and robustness against uncertainties p2| . We will prove 
that the proposed control law not only can ensure local exponential stability but also can give local finite-time stability. The 
exponential or finite-time stability can be switched by simply tuning a parameter in the control law. 

This paper is organized as follows. Section [II] formulates the formation control problem and presents our control law. We 



prove some useful lemmas in Section III The exponential and finite-time stability by the proposed control law is proved in 



Section IV Formation behaviors including collision avoidance between agents are analyzed in SectionjV] Numerical simulations 



are presented in Section |VT] A summary is given in Section |VII| 

II. Problem Formulation 

First of all, we present some notations that will be used through out the paper Let 1 = [1, . . . , 1]^ € K", and / be the 
identity matrix with appropriate dimensions. Denote | • | as the absolute value of a real number, and Null (•) as the null space 
of a matrix. Given x = [xi, . . . , x„]'^ G M" and p > 1, the p-norm of x is denoted as ||a; lip = (Er=i kd'')'^''- For the sake 
of simplicity, we omit the subscript when p = 2, i.e., denoting || • || as the 2-norm. Let be the entry at the ith row and jth 
column of a matrix, and be the zth entry of a vector. For a symmetric positive (semi-) definite matrix A, its positive real 
eigenvalues are denoted as < Ai(A) < \2{A) < ■ ■ ■ < Xn{A). 

Consider n (n > 3) agents in M^. The information topology of the agents is described by a graph Q = {V,£), which 
consists of a vertex set V — {1, . . . ,n} and an edge set £ C V x V. Each vertex in Q corresponds to an agent. If G £, 
then i and j are called to be adjacent, and i receives information from j. The set of neighbors of vertex i is denoted as 
Afi ^ {j E V : E £}■ A graph is undirected if each E £ implies E £, otherwise it is directed. A path from i 

to j in a graph is a sequence of distinct nodes starting with i and ending with j such that consecutive vertices are adjacent. 
If there is a path between any two vertices of a graph Q, then Q is said to be connected. An undirected circular graph is a 
connected graph where every vertex has exactly two neighbors. 
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Fig. 1: Illustration of a circular formation. 

An incidence matrix of a directed graph is a matrix E with rows indexed by edges and columns indexed by vertices 
Suppose {i, k) is the ith edge. Then the entry of E in the ith row and kth column is 1, the one in the ith row and jth column 
is —1, and the others in the ith row are zero. Thus by definition. El — 0. If the graph is connected, the corresponding E has 



rank n-1 |23 Theorem 8.3.1]. Then Null (E) = span{l}. 



A. Angle -constrained Target Formation 

In this paper we consider the distributed control of circular (or polygon) formations. The underlying information flow among 
the agents is described by a fixed undirected circular graph. By indexing the agents properly, we can have Afi — {i + 1, i — 1} 
for i £ {1, • • ■ ,n}, which means agent i can measure the bearings of agents i — 1 and i + 1. The indices i + 1 and z — 1 
are taken modulo n. Denote the position of sensor i as e M^, and the edge between agents i and i + 1 as = z^+i — Zi. 
The unit-length vector gi — ei/\\ei\\ characterizes the relative bearing between agents i + 1 and i (see Fig. [TJl. Hence the 
measurements of agent i consist of and — g^-i. It should be noted that agent i may measure gi and — g^-i in its local 
coordinate frame. While analyzing the dynamics of the entire system, we need to write these bearing measurements in a global 
coordinate frame. 

The angle subtended by agents i + 1 and i — 1 at agent i is denoted as 9i S [0,2tt). More specifically, rotating —gi-i 
counterclockwise through an angle di about agent i yields gi (see Fig.[T]i. Hence 9i and di+i are on the same side of edge 
for all is {1, . . . Consequently the quantity J2"=i invariant to the positions of the agents because the sum of the 
interior angles of a polygon is constant. In the target formation, the angle di is specified as 9* E [0, 27r). The target angles 
{d*}^^i should be feasible such that there exist {zi}^^^ (zi ^ Zj for i ^ j) to realize the target formation. Because the target 
formation is only constrained by angles, the realization will be non-unique. This paper makes no assumptions about parallel 
rigidity ||7|, p7) , p8) of the formation. As shown in Fig.|2] the angle-constrained target formation is invariant to the following 
motion patterns: rigid body transformation, scaling and edge parallel motion. The term invariant as used here means that these 
motions will not change the angles in the formation. Moreover, it is notable that = X]"=i regardless the motion of 

agents. 

'in some literature such as j23|, the rows of an incidence matrix are indexed by vertices and the columns are indexed by edges. 



4 




(a) Original fomiation. (b) After rigid body transformation. 




(c) After scaling. (d) After edge pai'allel motion. 



Fig. 2: The angle-constrained target formation is invariant to four motion patterns, (a) Original formation; (b) after rigid body transformation; (c) after scaling; 
(d) after edge parallel motion. 



B. Proposed Control Law 

The control task can be summarized as this: steer agents from their initial positions {zi{0)}f^i to final positions such that 
di — 9* for all i G {1, . . . , n}. For a proper definition of 9i, we need to assume Zi(0) ^ Zj{0) for all i ^ j. The feedback 
angle error of agent i is defined as 

Ei = COS 9i — COS 9* (1) 

= -gfgi-i - cos 6**. 

The reason why we use cosine functions to define the angle error Si is that cos 9i can be conveniently expressed as the inner 
product of the two bearing measurements gi and —gi-i- Suppose the motion model of each agent as a single integrator. The 
proposed control law for agent i is 

= sgn(ej)|ei|''(5r, - g^_i), (2) 
where a G (0, 1] and sgn(-) is the sign function defined by 

1 if > 

sgn(e,) = <( if = . (3) 
-1 ifei<0 

In the special case when a = 1, control law (|2]) becomes ij ~ £i{gi — ffi-i) because sgn(ei)|£j|° = e^. 

Remark 1: As will be shown later, control law (|2]i ensures local exponential stability if a = 1, and local finite-time stability 



if a <E (0, 1). Loosely speaking, finite-time stability means Si for all i converges to zero in finite time. See |22| or |24 Section 
4.6] for a formal definition of finite-time stability of nonlinear systems. In this paper, we will present a unified proof of the 
exponential and finite-time stability based on Lyapunov approaches. 

Remark 2: If a > 0, it is straightforward to see sgn(ei)|ei|'' is continuous in e^. If a = 0, control law (|2]l is discontinuous 
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in The discontinuous case with a = is out of the scope of this paper. One may refer to pS) , p6) for finite-time stability 
of discontinuous dynamic systems. 

Our proposed control law is inspired by the one in [T|, where the control law steers agents moving along the bisectors of 
the angles in a triangle. The velocity direction in control law (|2| is also along the bisector of Oi. But the feedback angle error 
and velocity are differently defined in (|2]i compared to |3|. It is noticed that gi — gi-i will vanish when 9i = tt. Hence the 
control law is ineffective in the case of Oi — tt even though is still nonzero. Moreover, when 9i — 0, agents i — 1 and i + 1 
are on the same side of agent i. Since bearing information is usually measured by optical sensors such as cameras, the bearing 
of agent i — 1 or « + 1 may not measurable by agent i due to line-of-sight occlusion in the case of 6i — 0. Therefore, we make 
the following assumption. 

Assumption 1: In the target formation, 0* ^ and 9* ^ vr for all i E {1, . . . ,n}. 

By Assumption [T] the angle 9* is in either (0,7r) or (tt, 27r). In other words, no three consecutive agents in the target 
formation are collinear. The collinear case is a difficulty in many formation control problems (see |[6), pTj , p3| , | [T6| for 
example). Because the angle error Si is defined using cosine functions and gi — gi-i vanishes in the case of 9i = tt, control law 
([2| inherently is not able to ensure global convergence to arbitrary feasible target formations. However, the benefit of control 
law (|2| is that its dynamics can be conveniently analyzed based on Lyapunov approaches. At last, we would like to emphasize 
that although gi and gi_i in (|2]i are expressed in a global coordinate frame, the control law can be implemented using agent 
i's bearing-only measurements in its local coordinate frame. 

III. Preliminary Results 

In order to analyze the stability by the proposed control law, we need to prove and introduce the following results. 
Lemma 1: Let A E M"^" be a positive semi-definite matrix with Xi{A) — and X2{A) > 0. An eigenvector with the zero 
eigenvalue is 1 = [1, . . . , 1]^ G M". Let 

U — {x E M" : ||a;|| — 1 and nonzero entries of x are not with the same sign}. (4) 

Then 

■mix^Ax^^. (5) 

x£U n 

Proof: By orthogonally projecting a: to 1 and the orthogonal complement of 1, we can decompose x as 

x = coXo + ciXi, (6) 
where ||a;o|| — \\xi\\ — 1, xq \\ 1, xi _L 1, cq = XqX and ci = xfx. Because Al — 0, I'^A — and xq || 1, we have 

x'^Ax = clxfAxi > clX2{A), (7) 
where the last inequality uses the well-known fact 

min xTAxi = X2(A). (8) 

xi_Ll,||a:i|| = l 

Denote Lp E [0, tt] as the angle between 1 and x. Then ci = sin (p. Geometrically speaking, the vector 1 is isolated from 
X E U hy the hyper-planes [x]i — with i E {1, . . . ,n}. Let pi E M" be a vector with zero as the zth entry and one as the 
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Others. Then is the orthogonal projection of 1 onto the hyper-plane [x]i = 0. Note the set U is open. The infimum angle 
(fiinf between 1 and x is the angle between 1 and pi. Thus I'^Pi ~ ||l||||pi|| cos(^inf, which implies cos(piiif = \/n — l/\/n 
and sin(^i„f = 1/y/n. The supremum angle i^sup between 1 and x is the angle between 1 and ~pi. In that case, we have 
costpsup = —\fn^-^l \/n and siniySgup = Therefore, 

inf{x : c\~ smip] — (9) 

substituting which into (|7]i yields ■ 
Lemma 2: Let x{t) be a real positive scalar variable for all t e [0, +oo). Given k E (0, 1), if the time derivative of x{t) 
satisfies 

\x{t)\ < exp ^ '^'^^) ' ^ ^ t'^' +'^)' ^^^^ 

then for all t e [0, +oo) has a finite upper bound. 

Proof: We will find an explicit finite upper bound by repeatedly calculating the integral in ( [TO] i and using the fact 

xit)<xiO) + J^\x{r)\dr. 

First of all, because a; > 0, we have —k/x < and hence 

< exp(O) = 1. (11) 

Then 

x{t)<x{0)+f ldT<t + c, c>x{Q). (12) 

^0 

Substituting the above inequality back into ( [T0| gives 

k 



\x{t)\ < exp f / ^dr 



fc 

c 



t + c 



(13) 



Then 



<Y^(^ + c)^-^ (14) 
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where the last inequality uses the fact c > a;(0), 1 — A; < 1 and hence a::(0) — c/(l — /c) < 0. Again substituting the above 
inequality into ( fTO] ) gives 



\x{t)\ < exp (^J^ 



exp(-l-^((t + cf-c'=)) 



ei-^exp(-^(t + c)M. (15) 



Denote = (1 - fc)/c'=. Then for all t e [0, +oo) 

x{t) < x{0) + e^-'' I e-^t^+^'^d 



T 







+ 00 



<x{0) + e^''' / e-^("+^) dr. (16) 



Let s = ii{t + c)*^ and hence dr = {l/k)ii i/'^gi/'^ ids. Then above integral becomes 

f+OO -1 /'+00 



/ e-^(-+^)'"dr= 

Jo ^ Jnc'' 



1 



flC 

4-00 



< ^A*"* / e-^s^J-Ms. (17) 







For any positive real constant a, the integral T{a) = J^°° e '^s" ^ds in fact is the Gamma function of a and has a finite 
value. Therefore, we conclude that x{t) for all t e [0, +oo) has a finite upper bound 

x{t) <x{0) + yfi'h^-''T (I) . (18) 



Lemma 3 ( /|2^ Lemma 2]): Let xi, . . . , a;„ > 0. Given p e (0, 1], then 



/n\P n /"\^ 
\i=l / i=l \i=l ) 

Lemma 4 ( / [2^ Corollary 5.4.5]): Let || • ||a and || • ||^ be any two vector norms on M". Then there exist finite positive 
constants Cm and C]\f such that C,„||a:;||Q < < CM||a;||a for all x £ M". 

IV. Stability Analysis 

In this section, we present a unified way to prove that the target formation by control law (|2]l is locally exponentially stable 
if a = 1, and locally finite-time stable if a e (0, 1). 

A. Basic Stability Analysis 

Denote e = [ei, . . . , £„]"^ G M" and z = [zf , . . . , zj^]^ G M^". From Q, it is straightforward to see e = implies i = 
and then e = 0. Hence e = is an equilibrium of the e-dynamics. We here omit giving the dynamic equation of e, whose 
dynamics, however, will be directly used in the following analysis. Now consider the Lyapunov candidate 



1 " 

i—l 
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Then V can also be written as V = l/(a + l)|l£|la+i where || • ||a+i denotes the (a + l)-norm. Clearly V is positive definite 
with respect to e = 0. In the case of a = 1, the Lyapunov candidate becomes V ~ l/2e^e, which is an ordinary quadratic 
function of e. 

If > 0, = ^ = (a + m = (a + l)sgn(e,)|£.r and lim,,_o+ - 0; if < 0, = ^^^^ = 

— (a + 1)(— Ei)'' — (a + l)sgn(ei)|ei|° and linig._>.o- — = 0. Therefore, for all Si we have 

^^^^-(a + l)sgn(£,)k.r. (21) 
For the sake of simplicity, denote ai — sgn(£i)|ei|° and cr = [cti, . . . ,(Tn]"^ G K". The time derivative of V is 



-E 



a + 1 

1 



1=1 



1=1 

n n 

1=1 i=l 

n n 

1=1 i=l 
n 

- ^(cTi+iffi+i + (Tig^^ifiji. (22) 



i=l 



Because t^i = ei/||ei||, we have 



9i ^-^-PjCj, (23) 
e,; 



where Pi = I — gigf . Note Pi is an orthogonal projection matrix satisfying Pf — Pi and P^ = P;. Moreover, Pi is positive 
semi-definite and Null(Pj) = span{(7i}. From control law (|2]i, we have 

= cri(.gi - (24) 

ii+i = o-i+i(gj+i - gi), (25) 

and hence 

= 0-^+154+1 + crj5i_i - (cTi+i + cri)5i. (26) 

Because P^g^ = 0, substituting the above into ( |23| ) gives 

1 



1 



PiicTi+igi+i + cngi^i). (27) 
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Substituting the above iji into ( |22| l yields 

" 1 

t=i ll^'ll 

Now we can claim the equilibrium £ = is at least Lyapunov stable. 



(28) 



B. Exponential and Finite-time Stability Analysis 

In order to further prove the exponential and finite-time stabiUty, we need the following results. 
For an arbitrary angle a, the rotation matrix 



R{a) 



cos a — sm a 
sin a cos a 



p2x2 



(29) 



rotates a vector in counterclockwise through an angle a about the origin. Thus it is clear that for all nonzero x G M^, 
R{a)x > when a E (— 7r/2,7r/2) (mod 27r); x^R{a)x = when a = ±tt/2 (mod 2tt); and x^ R{a)x < when 
a e (7r/2,37r/2) (mod 27r). Moreover, R^^{a) — R^{a) — R{—a) and R{ai)R{a2) = R{ai + 02) for arbitrary angles ai 
and a2- 

Lemma 5: Let g,^ = i?(7r/2)gi. It is obvious that ||.g,^|| = 1 and [gi^)'^ gi = 0. Furthermore, 

(i) p. = <?,^(5.^r. 

(ii) For i ^ j, {giYgj ^ -{gfYg^■ 

(iii) (3,^)^5^-1 > if e, e {0,tt); and {g^fg^-l < if 0, e (7r,2^). 

Proof: (i) Denote Gi = G M^^^. Then Gi is an orthogonal matrix satisfying GjGi — I. Hence we have 



g^gI + gt{gh^ ^G.Gf ^l. 

Thus gHgh"^ ^ I ~ g^gT = P^■ 

(ii) {gtVg, = fff i?^(7r/2)5, = glR{-7r/2)g, = gf R{-n)R{7r/2)g, = fff i?(-vr)5,^ - fff (-/).9,^ = 

(iii) By the definition of 9i, we have gi — R[9i){—gi-i) and hence gi^i — —R{—6i)gi. Then 



(30) 



Then it is straightforward to see the result in (iii). 

Because Pi — g^{g^)'^ as shown in Lemma |5](i), V in (|28]i becomes 



g, = - cos - 



= smf 



(31) 



V 



< 



"1 2 

- Juli ((5i^)^('^»+i.9i+i + <^igi-i)) 

1 " 

- v=^n — n— f7 X] ('^«+i(.9i^)^.9»+i + '7i(5./")'^5«-i)' 

Z^i=l II Si II .^^ 



cr2{giVg2 + criigiVgn 
f^iignVgi +<^n{gnVgn-i 



(32) 
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Because 



' (giVgn 












O'l 





(52^)^91 


(52^)^53 ... 







0-2 








{giVg2 ... 







0-3 










(gnVgn-i _ 







1 




-1 





-1 
1 



... 

... 

... 

1 



where the last equality uses the fact that {g^)^ 

V < 



{gifgn 






















(.92^)^51 












0-2 










(53^)^52 





















{g^Yc 


hi- 1 




















-1 = -(. 


li-ifgt as 


shown in 


Lemma [s] (ii), 


([32]i can be rewritten 



(33) 



-a^D'^E'^EDa. 



(34) 



We now are ready to present the main results of this paper 

Theorem 1: Under Assumption [T] the equihbrium e = is locally exponentially stable by control law (|2]i if a = 1, and 
locally finite-time stable if a e (0, 1). 

Proof: The proof consists of two steps. Step 1: prove D'^ E^ EDa is lower-bounded by KV°^ with iiT as a positive 
constant. Step 2; prove l/J27=i ll^^ll bounded from below by a positive constant. 

Step 1: 

Suppose £ ^ CT 7^ 0. Write 



cT^D'^E^EDa = 



a^D^E^EDa a^D^Da, 



a^DTDa 



(35) 



First, denote il{c) ~ {e : V{e) < c} with c > as the level set of V{e). Since < 0, the level set fl{c) is positively invariant 
with respect to (|2]). At the equilibrium point e = (i.e., 9i — 0* for all i e {1, . . . , n}), we have [D]ii — {gj')^gi-i 7^ because 
9* 7^ or TT as stated in Assumption [T] Thus by continuity there exists a sufficiently small c such that [D]ii 7^ for every point 
in n{c). Then D = D'^ is positive definite and hence Xi{D'^ D) > for all e e n{c). Because V = l/(a + l)||e||°+i[, the set 
il{c) is compact given a sufficiently small c. Hence there exists a lower bound Xi{D^D) > such that Xi{D^D) > Xi{D^D) 
and consequently 



(36) 
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for all e E ri(c). In addition, since 2a/ {a + 1) G (0, 1], we have 



, 1 \ STT " 

^ — T El^^l'" (by Lemma|3j 



a + 1, . , 
^ i— 1 

1 \ " + 1 



a + 1 



la 

1 \ » + ! 



a + 1 



Y^a\ (by = 

=1 

cr^cr. (37) 



From ([36| and ( (37] l we have 



= (a + l)^Ai(i:>^D) (38) 

for all e e fi(c) \ {0}. Because Xi{D'^D) = Ai(i:i2) ^ min,((g,^)^5,_i)2 < 1, we have Ai(D^D) < 1. 

Second, by the definition of incidence matrices, E is an incidence matrix of a directed and connected circular graph. By p3] 
Theorem 8.3.1], we have rank(£^) = n — 1 and consequently ra.nk{E^ E) = n—1. Note El = 0. Then Null (E'^E) = spanjl}. 

Denote 5i = 9i — 9* and 

cos dj — cos 9* 

Hence Si = WiSi. Since lim^.-^g* Wi — —sm9*, the equation Si — wiSi is still valid even if 9i — 9* = 0. Then we have 

e = WS, (40) 

where W — diag{wi, . . . , w„} G M"^" and 6 = [Si, . . . , 6n]'^ G M". There exists sufficiently small c such that 9i{0) is 
sufficiently close to 9* and hence 0^, ^* e (0, tt) or 0^, 0* e (tt, 27r) for all e G fi(c). Clearly < when 9^, 9* £ (0, tt), and 
> when 6'^,^* e (tt, 27r). Recall [g^)'^ g-i-i > when e (0,7r), and {g^Y gi^i < when 9i £ (tt, 27r) as shown in 
Lemma |5](iii). Thus {g^)'^ gi^iWi < for all i G {1, . . . ,n} and consequently the diagonal entries of DW are with the same 
sign. Moreover, because 9i = 9*, the nonzero entries in 6 are not with the same sign. Therefore, the nonzero entries 
of De = DWS are not with the same sign. Furthermore, because ai have the same sign as e^, the nonzero entries of Da are 
not with the same sign either By Lemma [T] we have D<t/\\D<t\\ e U and 

Because J27=2 '^ii^'^ — ir{E'^E) = — 2"- and \2{E'^ E) is the smallest positive eigenvalue of E'^E, we 

have (n - 1)A2(-B^£') < 2n. Thus X2{E'^E)/n < 2/{n - 1) < 1 as n > 3 and the equality holds only if n = 3. 
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Substituting ([38]) and (|4T]) into ([35]) and (|34]) yields 

1 



V < — r—ra'^D'^E'^EDa 



Er=i 



<- ^»\| (42) 



where 



K = {a + l)—^\AD^D)-^ '-. (43) 

n 

As a e (0, 1], (a+ 1)^ < 2. Recall X^{D'^D) < 1 where the equality holds only if {g^Y gi^i = ±1, i.e., gi ± for all 
i e {1, . . . , n}; and \2{E'^ E)/n < 1 where the equality holds only if n = 3. Obviously gi _L for all i will not hold for a 
triangle with n = 3. Hence \2{E'^E)/n = 1 and \^{D^ D) = 1 cannot hold simultaneously. Then X^{D'^ D)X2{E'^ E) /n < 1 
and hence K e (0, 2). 
5fep 2.- 

Since the distances between agents are not controlled directly, we cannot simply rule out the possibility that J27=i W'^iW 
may go to infinity. Next we will prove J27=i ll^^ll bounded from upper by a finite positive value and hence l/J27=i 11^*11 
bounded from below by a positive constant. 

Recall Ci ~ o'j+i(.9i+i — gi) + cnigi-i — gi) as shown in ( p6] ). Denote p — X]"=i ll'^ill- Then the time derivative of p is 

d||e,|| 



E 

i=l 



dt 



i=l 
n 

= ^9l[<^t+i{9t+i - ft) + cri{gi-i ~ ft)] 

i=l 
n 

= Et^'+if^^f ft+1 - 1) + cr«(ft^ft-i - 1)] 

2 = 1 

n n 

= E '^''+i(ft^ft+i ~ 1) + E ^^idJgi-i - 1) 

i=l j=l 

n n 

= E '^^ial-igi - 1) + E ^»(ft^ft-i - 1) 

i=l i=l 

n 

= 2^a,(5fg,_i-l) 

1=1 

= 2D'^cr, (44) 
where v — [vi, . . . , Vn]'^ £ M" with Vi = gfgi~i — 1. Then by the Cauchy-Schwarz inequality 

|p| = 2|«^a|<2||i;||||a||<2/3||a||, (45) 

where /3 is the maximum of ||u|| over the compact set f2(c). 

For any c < 1, we have < 1 and then > F as 2a/(l + a) < 1 for all e G ^l{c). Thus 

K 2a K 

V < < V. (46) 

P P 
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By the comparison lemma pS\ Lemma 3.4], the above equation implies 



* K 



Vit) < V{0) + I --^V{T)dT. (47) 







Applying the Gronwall-Bellman inequality |29 Lemma A.l] to the above equation gives 



In addition, we have 



Vit) <ViO)e,,p{ 1^ -^dr). (48) 



/ 1 \ STi 1 " 

^ — T ^ E l^^l"" Lemmas 



i=l 

2a 

1 \ »+l 1 



a + 1 



<rf, (49) 



which implies 



Substituting (|50]) and (|48]) into (|45]) yields 



< (a + (50) 



IpI < 2/3|k|| 

< 2/3VKt/(0)* exp ^ -^i^dr^ 



< exp ( / -^ii^dr ) . (51) 



—K 
P(^) 

The last inequality uses the fact that there exists a sufficiently small c such that 2f3y/KV{0)^+^ < 1 when V{0) < c. Because 
a/(l + a) e (0, 0.5] and K e (0, 2), we have aK/{a + 1) G (0, 1). Then by Lemma |2j /? is bounded from upper by a finite 
constant. Denote this upper bound as 7. Then V in ( |42| i becomes 

K 2a 

V< V^. (52) 

7 

For ease of presentation, we make a temporary assumption here that no collision between any agents occurs for all t £ 
[0, +00). Without this assumption, inequality ( |52| i only holds until collision happens. In the next section we will prove that 
the proposed control law ensures collision avoidance. Then this assumption can be removed. Under the collision avoidance 
assumption, inequality ( |52] l holds for all t E [0, +00). 



If a € (0, 1), 2a/ (1 + a) G (0, 1). By |22 Theorem 4.2], the solution to (|2]i starting from f2(c) converges to the equilibrium 
e = in finite time. 



If a = 1, 2a/(l + a) = 1. By |24 Theorem 3.1], the equilibrium e = is locally exponentially stable. 
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V. Formation Behaviors in the Plane 

Because the target formation is only constrained by angles, the positions of the agents and the inter-agent distances are not 
controlled in the final converged formation. In addition to the dynamics of e as analyzed in the last section, it is also important 
to investigate the behavior of the positions of the agents in the plane. 

From control law (|2|, one trivial behavior of the formation is that i = if £ = 0. In other words, we can rule out the 
possibility that the agents are still moving while the target angles have been achieved. Another trivial behavior of the formation 
is that i = if e 0. That is because = if and only if e 0. Intuitively speaking, it is impossible that control law (|2]) 
only changes the orientation, translation or scale of the formation while keeping all angles unchanged. 

Collision avoidance is an important issue in all kinds of formation control. It is especially important for formation control 
using bearing-only measurements because the distance between any two agents cannot be measured or controlled directly. 
The next result shows a simple but important behavior of the formation by the proposed control law: when the initial angles 
are sufficiently close to the target angles, the final converged positions of the agents are also sufficiently close to the initial 
positions. As a consequence, the proposed control law can locally ensure collision avoidance. Then the temporary assumption 
on collision avoidance in the proof of Theorem [T] can be removed. 

Theorem 2: Suppose in the initial formation Zj(0) 7^ Zk{Q) for all j, fc e {!,..., n} and j 7^ k. Then there exists a 
sufficiently small ||£(0)||a+i and a positive constant 77 such that the distance between z(t) and z(0) as shown below satisfies 

n 

^||z,(t)-z,(0)|| <ry||e(0)|U+i (53) 
1=1 

for all t E [0, +00). As a result, there exists a sufficiently small ||£(0)||a+i such that collision avoidance between any agents 
can be ensured by control law (|2]). 

Proof: We prove by contradiction. Suppose no collision between any agents occurs until T* E (0, +00). As shown in the 
proof of Theorem [T] there exists a sufficiently small c > such that inequality ( |52] i holds for t E [0, T] with T < T* if £(0) E 
n{c). When c is sufficiently small, we also have V{0) < 1 and hence V{t) < 1 for all t E [0,T]. Since 2a/(l + a) E (0, 1], 
we have >V and consequently 

K 2a K 

V < < V, (54) 

7 7 

which implies 

V{t) < V{0)e~^* (55) 
for all t E [0,T]. Substituting V = l/(a + 1)||£||°|} into (|55]) yields 

\\e{t)\\a+i < ||£(0)|U+ie-(^*. (56) 
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Thus the following "distance" between z{t) and z(0) satisfies 



;t II, 1,1 

^||z,(t)-z,(0)|| / £,(5,-5,_i)dT 

<E / |£.|ll5^-5.-i||dr 

.-1 "'0 



(by control law ([2|) 



1=1 

/•t n 

< 2 / V |e,|dr (becuause - < \\g,\\ + \\g,_i\\ = 2) 

•^0 ^=l 

= 2 f \\e\\idT 
ft 







<2C [ ||e||a+idr (by Lemma |4| 
Jo 



<2Cj ||e(0)||„.+ie-(™"dr (by Q) 
= 2C\m\\aJ-^{l-e-r^^) 

<?^^lk(0)IU.i (57) 




for all t e [0,T]. 

Now suppose agents j and k collide at time T* . Then Zj{T*) — Zk{T*) and 

n 

J2 P^^iT*) - ^.(0)11 > \\z,(T*) - z,{0)\\ + \\z,{T*) - Zfe(0)|| 

i=l 

> ||z,(r*)-z,(0)-Zfc(T*) + Zfc(0)|| 

= ||zfc(0)-z,(0)||. (58) 

However, given any positive constant r with r < |lzfe(0) — Zj(0)||, by ( |57| ) we can always find a sufficiently small e(0) such 
that 

n 

^||z,(i)-z,(0)|| <?7||e(0)|U+i 

i=l 

< ||zfc(0)-z,(0)|l -r (59) 



for all t e [0,T] and all T < T*. Clearly ( |58| l conflicts with ( |59l ). Therefore, collision avoidance can be guaranteed by control 
law (|2]i if e(0) is sufficiently small. Consequently ( |52| ) and ( |57| i holds for all t e [0, +oo). ■ 

VI. Simulations 

In this section, we present numerical simulations to illustrate our theoretical analysis. 

The target formation in Fig. [3] is a non-convex star polygon with n = 5. The angle at each vertex in the target formation 
is 9* ~ 36 deg for all i £ {1, . . . , 5}. From an initial formation, the proposed control law effectively reduces the angle errors 
and steers the formation to a target formation. Because the target formation is constrained only by angles, the shape of the 
final formation does not look like a regular star polygon though all angle errors achieve zero. 

The target formation in Fig. |4]is a ten-side polygon, where the angle at each vertex is 6* — (10 — 2) x 180/10 = 144 deg 
for all i e {1, . . . , 10}. In the stability analysis, we have mentioned there exists a sufficiently small c such that 9i,0* E (0, tt) 
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a=1 




(a) Formation evolution witli a = 1. 




(c) Formation evolution with a = 0.3. 




40 60 
Time (sec) 

(b) Angle en'or and Lyapunov function evolution with a = 1. 
8=0.3 



100 




40 60 
Time (sec) 



100 



(d) Angle error and Lyapunov function evolution with a = 0.3. 



Fig. 3: Control results by the proposed control law with n = 5 and 9'^ 



: e* = 36 deg. 



or (tt, 27r) for all points in rt{c). However, it is notable that in the initial formation 6',;(0) = tt for i = 2,3,7,8,10 and 
05(0) G (7r,27r) but 9'^ E (0,7r). The target formation can still be achieved. Hence from the simulation, the attractive region 
of the target formation by the proposed control law is not necessarily small. 

Figures |3] and |4] show that the angle errors and the Lyapunov function converge to zero in finite time if a < 1. 

VII. Conclusion 

This paper addresses a relatively new formation control problem: distributed control of angle-constrained circular formations 
using local bearing measurements. The proposed control law ensures exponential or finite-time convergence. The presented 
Lyapunov analysis might be helpful for future research on more complicated target formations. 

The work in this paper is the first step towards the control of generic angle-constrained formations using bearing-only 
measurements. Many interesting problems in this field are still unsolved. In the future, one immediate research plan is to study 
the control of generic angle-constrained target formations, in which one agent may correspond to multiple constrained angles. 
It is also meaningful to study formation control involving agent position estimation based on bearing measurements. 
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a=1 




(c) Formation evolution with a = 0.6. 
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(b) Angle en'or and Lyapunov function evolution with a = 1. 
8=0.6 
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(d) Angle error and Lyapunov function evolution with a = 0.6. 



Fig. 4: Control results by the proposed control law with n = 10 and SJ' 



: 6»* = 144 deg. 
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